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Update of MRST partons 
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We discuss topical issues concerning the MRST parton analysis. We introduce some new parton sets in order 
to estimate the uncertainty in the cross sections predicted for W and Z production at the Tevatron and the LHC. 
We compare with CTEQ5 partons. 



MRST is a standard next-todeading order 
global analysis [ Q of deep inelastic and related 
hard scattering data. The partons were presented 
at DIS98 and the scheme dependence and higher 
twist effects studied in a subsequent analysis [ || . 
Here we concentrate on topical issues. 

1. Light quark distributions u,d,s 

u and d are pinned down by DIS data, with the 
slope of d/u also being constrained for x < 0.3 
by the CDF W^ rapidity asymmetry data. In 
fact d/u is determined up to i ~ 0.7 by the 
NMC F£ / F% data, provided that we assume that 
there are no deuterium binding corrections for 
x > 0.3. Yang and Bodek [ || included such a 
binding correction and described the data with 
modified partons which satisfy d/u — ► 0.2 as 
x — * 1, whereas MRS(T) and CTEQ parametriza- 
tions have d/u — ► 0. We repeated the global anal- 
ysis, forcing d/u — > 0.2 and found an equally good 
fit, apart from the description of NMC Fg/F% for 
x > 0.5. To fit these latter (uncorrected) data we 
require d/u — > 0, see also ref. [ ||. We note, that 
as statistics improve, the HERA charge-current 
DIS data will directly constrain d/u at large x. 
So far the preliminary data and the MRST pre- 
diction are in agreement. 

d/u is determined for 0.05 < x < 0.3 by the 
E866 pp,pn — » u + \i~ X data, and also by the 
HERMES data for ep, en -> e^X. The MRST 
partons are in agreement with both the new E866 
[ | and HERMES [ § data. 



In principle v and v DIS data can determine 
s and 5 separately. However we use the CCFR 
dimuon data to determine s + s. To quantify the 
uncertainty we obtained two new sets of MRST 
partons, s f and s J, corresponding to (s + s)/(d + 
u) = 0.5 ± 0.06 at Q 2 = 1 GcV 2 , which embrace 
the spread of these data. 

2. Heavy quark distributions c, b 

For a heavy quark, H = c (or b) we should 
match, at Q 2 = m 2 H , the fixed flavour num- 
ber scheme, describing jg — > HH with rif = 
3(4), to the variable flavour number scheme with 
rif + 1 flavours for Q 2 > m 2 H . MRST use the 
Thorne-Roberts prescription [ pi] in which F^ 
and dF^ 1 /d\nQ 2 are required to be continuous 
at Q 2 = m 2 H . (CTEQ use an alternative, ACOT, 
prescription.) In this way the charm distribu- 
tion is generated perturbatively with the mass of 
the charm quark m c as the only free parameter. 
Not surprisingly, due to the g — > cc transition, 
the charm distribution mirrors the form of the 
gluon distribution. Again, we have obtained two 
new MRST sets, c f and c J,, corresponding to 
m r = 1.35 ±0.15 GeV. 



3. The gluon distribution g 

In principle many processes are sensitive to 
the gluon distribution, but it is still difficult to 
nail down for x ^ 0.2 where it becomes increas- 
ingly small. In this region MRST used the WA70 



prompt photon data to determine the gluon. 
However this process suffers from scale depen- 
dence (reflecting the higher order corrections), ef- 
fects of intrinsic kx (needed to describe the E706 
prompt photon data) and uncertainties due to 
fragmentation effects and isolation criteria. To 
explore an acceptable spread of gluon distribu- 
tions at large x, MRST presented three sets (g |, 
default MRST set, g |) corresponding to aver- 
age intrinsic kr — 0, 0.4 and 0.64 GeV respec- 
tively for the WA70 data. The arrows indicate 
the relative sizes of the gluons in the large x re- 
gion (x > 0.2). The corresponding (kx) needed 
to obtain a good description of the E706 prompt 
photon data is about 1 GeV. 

Due to the uncertainties associated with 
prompt photon production, the recent CTEQ5 
analysis [ || omits these data and instead deter- 
mines the large x behaviour of the gluon using 
the single jet inclusive Et distributions measured 
at the Tevatron. The above three MRST sets of 
partons describe the shape of the jet Et distri- 
butions, but need to be renormalized upwards by 
7%, 13%, 17% respectively to give the correct nor- 
malization, with g | being closest to CTEQ5. 

Turning now to the small x domain, the gluon 
is well constrained by the observed behaviour of 
dF 2 /d\ogQ 2 . To be precise the HI and ZEUS 
data determine the gluon for x < 0.01, and the 
NMC data for x < 0.1. Inspection of F% at low 
x, low Q 2 shows that the data have flattened 
out, indicating a smaller value of dF 2 jd log Q 2 ~ 
c t s(Q 2 )Pqg®g- Also as increases with decreasing 
Q 2 , and as a consequence g becomes valence-like 
for Q 2 ~ 1 GeV 2 , which is about the lowest value 
of Q 2 for a satisfactory DGLAP description. 

4. W production at the Tevatron and LHC 

The total cross sections for W and Z hadropro- 
duction are known to NNLO [ [| and the input 
electroweak parameters are known to high ac- 
curacy. The main uncertainty comes from the 
input parton distributions and, to a lesser ex- 
tent, ag [ [Hj. For example, the results for W 
production at the Tevatron from the above sets 
of partons are shown in Fig. 1. In addition we 
show the cross sections obtained from parton sets 
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Figure 1. NNLO [ |9| predictions of the W pro- 
duction cross section at the Tevatron from seven 
different sets of NLO MRST partons, and from 
CTEQ5HQ partons. Also shown are the cross 
sections measured by the CDF and DO collabora- 
tions. 



a s Tti a s II corresponding to global analyses 
with as(Aff) = 0.1175 ± 0.05, and from parton 
sets q f, q J, corresponding to global fits with the 
HERA data renormalised up, down by 2.5%. The 
differences coming from s "f, s J, and from c f, 
c I are very small (and so the predictions are 
not shown in Fig. 1), since in the evolution up to 
Q 2 = Myr the transitions g — ► ss, cc wash out the 
dependence on the s, c input. 

We see that the values from all the MRST sets 
are within about ±2.5% of the MRST default 
set. The bounds are actually given by the q f 
and q I sets. At first sight we might anticipate 
these bounds would be ±5% since qq initiates W 
production. However the difference due to using 
q t and q J, is diluted by the evolution from Q 2 
(HERA data) -> Q 2 = M 2 ,. 

We also show in Fig. 1 the prediction obtained 
using CTEQ5HQ partons. Surprisingly it differs 
by about 8% from that of MRST. Since we hope 
that the rate of W production is accurately pre- 
dicted at the Tevatron and at the LHC, it is cru- 
cial to understand the origin of this difference. 



5. Comparison of MRST and CTEQ5 

Differences between the MRST and CTEQ par- 
tons can arise from 

• the data sets included in the global fit, the 
Q 2 , W 2 . . . cuts used, 

• the normalization of the data sets, 

• the treatment of heavy flavors (massless, 
TR and ACOT schemes), 

• differences in the evolution codes. 

Inspection of the resulting partons show that 
CTEQ5 has, first, a larger gluon at large x than 
MRST (due to the different data fitted), second, 
a larger charm distribution, and finally a smaller 
d/u at very large x. 

The unexpectedly large difference between the 
values of the W cross section calculated from the 
MRST and CTEQ5 partons, seen in Fig. 1, ap- 
pears to be due to the more rapid evolution of 
the sea quarks of CTEQ5. A similar effect was 
found in a comparison of various NLO evolu- 
tion codes made during the 1996 HERA work- 
shop [|llj. Neither CTEQ, nor MRS, agreed com- 
pletely with the standard code after evolution up 
to Q 2 = 100 GeV 2 . The CTEQ sea was some 4% 
high and the MRS sea was 1% low and the gluon 
about 1.5% low. A careful check of the MRS code 
has recently revealed a small error in the NLO 
contribution to the evolution of the gluon. After 
this correction the agreement between MRS and 
the standard evolution code is exact. The cor- 
rected code increases the quarks in the relevant 
x and Q 2 range by typically 1%, and hence the 
MRST values of the W cross section at the Teva- 
tron shown in Fig. 1 are increased by about 2%, 
so that 

B w cr w = 2.45 ± 0.06 nb. (1) 

A full discussion will be given in ref . [ |l(J . 

6. Outstanding problems 

At present the two major deficiencies in the 
parton determinations are the absence of a com- 
plete NNLO analysis and the need to determine 



realistic errors on the partons. This last task is 
especially problematic. Determinations in certain 
(x, Q 2 ) domains often rely dominantly on one ex- 
periment. Moreover we need some way to quan- 
tify the uncertainties due to kx smearing, heavy 
target and other nuclear corrections, higher twist 
cuts etc. However, as evidenced by the contri- 
butions at this conference, continued progress is 
being made in improving our knowledge of par- 
tons. 
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